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Abstract 

In recent experiments of laser pulse interaction at relativistic intensities with a low density 
plasma, the proton radiography technique showed evidence of long-lived field structures generated 
after the self-channeling of the laser pulse. We present 2D particle- in-cell simulations of this 
interaction regime, where the dynamics of similar structures has been resolved with high temporal 
and spatial resolution. An axially symmetrical field pattern, resembling both soliton-like and 
vortex structures, has been observed. A study of the physics of such structures and a comparison 
with experimental data is in progress. 
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I. INTRODUCTION 



The study of the propagation of intense laser pulses in low density plasmas is relevant 
to several highly advanced app.ications, .ae.ud.ng electron and ion acce.eration de- 
velopment of X- and 7-ray sources, and fusion neutron production [3J, |4||. It is also of 
fundamental interest, due to the variety of relativistic and nonlinear phenomena which arise 
in the laser-plasma interaction [5j. 

During the interaction of high-intensity laser pulses with matter, strong transient electric 
fields result from the charge separation induced either by the ponderomotive force of the 
laser or by the instantaneous flow of hot electron currents. Depending on the interaction 
conditions, these fields can reach extremely large amplitudes, up to the TV/m range. Electric 
fields due to charge separation can drive the expansion of the ions of the plasmas, leading 
to production of multi-MeV ion beams in interaction with thin foils , or to Coulomb 
explosion of plasma channels in interaction with underdense plasmas. A major step forward 
in the detection of such fields has been marked by the development of the proton imaging 
and deflectometry techniques which, employing laser-driven protons as a particle probe, have 
proven to be an exceptionally useful tool for the investigation of ultra-fast plasma dynamics. 
The proton beams emitted from a laser - irradiated foil are highly laminar [§| and, for 
projection purposes,can be described as emitted from a virtual point-like source located in 
front of the target [9|. A point-like projection imaging scheme is therefore automatically 
achieved. The unique capability of this technique to detect electrostatic fields in plasmas 
has allowed the retrieval of direct information on electric fields arising through a number of 
laser-plasma interaction processes. The high temporal resolution, related to the picosecond 
duration of the proton burst at the source is fundamental in allowing the detection of highly 
transient fields following short pulse interaction. 

In this paper we present results from some recent experiments as well as the results of 
the simulations obtained for the range of numerical parameters close to the experiment. In 
these investigations, proton probes have been used to investigate ultrafast plasma dynamics 
initiated by high-intensity (10 18 -=-10 19 W/cm 2 ) interaction, by detecting the ensuing electric 
and magnetic fields, as well as to give evidence of as long-lived field structures generated 
after the self-channeling of the laser pulse. In particular, we discuss the case of laser pulse 
propagation through underdense plasmas. Beside providing novel information on the physics 
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of this interaction regime, the data provides a clear example of the diagnostic capabilities 
of probing techniques employing laser-driven multi-MeV protons. 



II. EXPERIMENTAL SETUP AND SIMULATION APPROACH 

The experiment was carried out at the Rutherford Appleton Laboratory, employing the 
VULCAN Nd-Glass laser system [lo|, providing two Chirped Pulse Amplified (CPA) pulses, 
with 1.054/im wavelength, synchronized with picosecond precision. Each of the beams de- 
livered approximately 30 J on target in 1.3 ps (FWHM) duration. By using off-axis parabo- 
las, the beams were focused to spots of lO/mi (FWHM) achieving peak intensities up to 
3xl0 19 W/cm 2 . The short pulses were preceded by an Amplified Spontaneous Emission 
(ASE) pedestal of 300 ps duration One of the beams was directed to propagate through 
He gas from a supersonic nozzle, having a 2 mm aperture, driven at 50 bar pressure. The 
interaction was transversely probed by the proton beam produced from the interaction of 
the second CPA beam with a flat foil (a lO/mi thick Au foil was typically used), under the 
point projection imaging scheme (see Fig. [1]). Even if the targets nominally do not contain 
hydrogen, hydrocarbon contaminants that are always present on the target surfaces ensure 
a suitable source of protons [l^ . The detector was a multilayered radiochromic film detec- 
tor placed at distance 2 3 cm from the gas jet. In the condition of the experiment, this 



provided a multi-frame temporal scan of the interaction for up to 50 ps in a single shot 13]. 
The temporal resolution of each frame was typically of order of few picoseconds. 

For better understanding of the processes observed in the laboratory experiment particle- 
in-cell (PIC) simulations were carried out. The numerical model is based on Vlasov equations 
for both electron and ions and Maxwell equations for the electromagnetic field. Schematic 
drawing of the the laser-plasma interaction geometry for two-dimensional (2D) PIC simu- 
lations is shown in Fig. [2j The reduction to a 2D geometry is necessary because a fully 
3D simulation with both proper numerical resolution and spatial and temporal scales equal 
to the experimental ones is far beyond present-day supercomputing power. Apart from 
this unavoidable limitation the simulation parameters are close to the experimental ones. 
The simulations were performed on a Linux cluster at the CINECA supercomputing facility 
(Bologna, Italy). 
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III. CHANNELING IN UNDERDENSE PLASMA 



The data presented in Fig. [3] refers to the early stage of the interaction and clearly show 
that a charge displacement channel is produced by the laser pulse. The laser pulse propagates 
from left to right. 

A white channel with dark boundaries and bullet shaped leading edge is clearly visible in 
the frames (a) and (b). The channel propagates along the axis with the velocity that is equal 
to the speed of light within the experimental error. In the trail of the channel the proton flux 
distribution around the axis changes qualitatively (see Fig. [3b), showing a dark line along 
the axis. The white channel indicates the presence of a positively charged region around the 
laser axis. In this region the electric field points outward in the radial direction. The central 
dark line (labeled as "III" in Fig. [3]) observed at later times in the channel suggests that at 
this time the radial electric field must change its sign at some radial position leading to the 
focusing of the probe towards the axis. 

2D electromagnetic PIC simulations in planar geometry qualitatively reproduce the ex- 
perimental results. In the simulations the laser pulse was s-polarized (in the experiment 
the interaction pulse was polarized along the axis of proton probing), thus E z represents the 
amplitude of the propagating electromagnetic pulse and E y is generated by the space-charge 
displacement . 

The simulation clearly shows the formation of the leading part of the laser pulse. This 
force generates the radial space-charge electric field, which points away from the laser prop- 
agation axis. In the channel trailing edge one observes two narrow ambipolar fronts, slowly 
moving away from the axis. As was shown in Ref. [14] such a radial electric field profile 
produces a pattern in the proton images similar to that observed in the experiment. 

The essential dynamics of the formation of the charged channel and the subsequent ion 
evolution could be reproduced also by using a simple electrostatic, ponderomotive model 
in ID cylindrical geometry The simulation of proton deflection in the field patterns 

reproduced by such model gave proton images very similar to the experimental ones, thus 
validating the approach for a study of ponderomotive ion acceleration in the channel. Fig lVIl 
gives snapshots of the electric field and density profiles and the ion and electron phase space 
from a typical ID simulation. An interesting effect is the "echo" behavior of the electric field, 
which almost vanishes at the end of the laser pulse, to reappear later at a point where the ion 
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fluid, accelerated during the laser pulses, forms a sharp spike and undergoes hydrodynamical 
breaking and lead to electron heating. The ambipolar field structure can be thus interpreted 
as the sheath field generated due to electron heating near the density spike. More details 
can be found in Ref. jl5| . 

IV. STRUCTURE FORMATION 

At later times, typically on the order of 6-8 ps, the development of quasiperiodic modu- 
lations inside the channel was observed in the experiment (Fig. ED- These structures evolved 
into circular structures which were observed to decay on hydrodynamic time scales. These 
structures are tentatively being interpreted as related to the growth of electromagnetic soli- 
tons inside the channel. " Soliton-like" structures form in high-intensity interactions due to 
trapping of the red-shifted electromagnetic radiation by the ambient plasma, which behaves 



as over-dense for them [l6j. The ponderomotive force associated to the trapped radiation 
expels the electrons from the core producing a positively charged sphere, which can deflect 
the protons of a probe beam. Solitons (or post-solitons, i.e. larger structures arising from 



soliton merging 17|) can then be imprinted over the RCF as a "white" region 181 ] . 

The simulations, al later times, give evidence of the formation in the lower density region 
of a deep channel (due to ponderomotive expulsion of electrons and ions out of the axis). 
In the higher density region, a variety of nonlinear effects such as laser beam breakup into 
filaments and soliton formation occurs (Fig. [7]). 

In the wake of the laser pulse (see Fig. [8]) we observe the formation of "standing" (low 
frequency and almost zero propagation velocity), structures on the two sides of the laser- 
drilled channel, distributed as two (anti) symmetrical rows with respect to the x-axis. Their 
generation and dynamics may be due to the excitation of surface modes localized on the 
channel walls. In fact, the dispersion relation of surface waves allow the propagation of low- 
frequency electromagnetic radiation from the interaction region down to the low-density 
plasma. Here, due to the decreasing density, the group velocity of surface waves slows down 
and the field structures become quasi-stationary. 

In the direction perpendicular to the simulation plane these structures appear to have 
both a steady magnetic field B z (so that they resemble a vortex row) and a steady electric 
field E z localized into a density depression (analogous to zero- frequency "solitons"). From 
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the evolution of the ion density and of the electrostatic field (E x ,E y ) in the simulation 
plane we infer that the late evolution of these structure is due to local ion acceleration long 
their "radius". Frequency analysis of fields and currents gives evidence of the generation of 
slowly- varying fields. 

V. CONCLUSION 

In conclusion, we have reported the direct experimental study of electric field dynamics in 
a charge displacement channel produced by the interaction of high intensity laser pulse with 
an underdense plasma. The use of laser-driven proton probes as a diagnostic tool has proven 
crucial in obtaining the information on the structure and development of these fields with 
picosecond resolution. The experimental investigations are supported by massively parallel 
PIC simulations and simplified analytical/numerical modeling able to addressing the "long" 
picosecond timescale. 

At long times, both experiment and simulation give evidence of long-lived, "quasi- 
periodic" field structures in the plasma channel. We observe both sub-cycle electromagnetic 
"solitons" at the laser pulse front, and quasi-periodic structures on the channel sides in the 
pulse wake. These latter structures are reminiscent of both vortex and solitonic structures. 
The spatio-temporal evolution and the frequency spectrum of these structures have been 
characterized in detail, providing a rich information for the theoretical interpretation. 
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FIG. 1: Experimental arrangement. 



FIG. 2: Initial configuration and the model density profile (left frame). The electron density 
profile (right frame) measured along the propagation axis before the high-intensity interaction is 
consistent with the neutral density profile of the gas jet and indicates complete ionization by the 



ASE prepulse 



14] 



FIG. 4: Distributions of ions density and E z and E v electric field components obtained from PIC 

n 

simulations at 2 ps [141 ] . The bottom plots show the lineout of the radial electrostatic field at two 
positions along the x axis. 



FIG. 5: Simulation results from a ID electrostatic, ponderomotive model. The profiles of elec- 
tric field E r (blue, thick line) and ion density (red, dash-dotted line), and the phase space 
distributions of ions fi(r,p r ) and electrons f e (r,p r ) are shown at various times. 



FIG. 6: Proton projection images of the interaction region taken at: (a) 9 ps and (b) 45 ps after 
the pulse propagation through the center of the frame |l9j]. 



FIG. 7: Distribution of the ion density obtained from PIC simulations at 2.83 ps. 



FIG. 8: Distribution of the ion density (top frames), B z component of the magnetic field (middle 
frames) and E z component of the electric field (bottom frames) obtained from PIC simulations at 
different times. 

FIG. 3: Proton projection images of the interaction region at different times jl4]. The time labels 
give the probing time of the proton propagating along the probe axis. White regions correspond 
to lower proton flux, meanwhile dark regions correspond to higher proton flux. 
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